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Abstract. We discuss contact invariant structures on the space of solutions of 
a third order ordinary differential equation. Associated to any third-order dif- 
ferential equation modulo contact transformations, Chern IChe40l introduced 
a degenerate conformal Lorcntzian metric on the space of 2-jets of functions 
of one variable. When the scalar invariant of Wiinschmann [WiinOSl vanishes, 
the degenerate metric descends to a proper conformal Lorentzian metric on 
the space of solutions. In the general case, when the Wiinschmann invariant is 
not zero, we define the notion of a causal geometry, and show that the space 
of solutions supports one. The Wiinschmann invariant is then related to the 
projective curvature of the indicatrix curve cut out by the causal geometry 
in the projective tangent space. When the Wiinschmann vanishes, the causal 
geometry is then precisely the sheaf of null geodesies of the Chern conformal 
structure. We then introduce a Lagrangian and associated Hamiltonian from 
which the degenerate conformal Lorentzian metric are constructed. Finally, 
necessary and sufficient conditions are given for a rank three degenerate con- 
formal Lorcntzian metric in four dimensions to correspond to a third-order 
differential equation. 

1. Introduction 

The purpose of the present paper is to study the geometry of third-order ordinary 
differential equations: equations of the form 

y'" ^F{x,y,y',y"). 

By setting p ~ y' and q = y" , the solutions of this equation coincide with curves in 
the second jet space with coordinates {x,y,p,q) that are everywhere tangent to 
the contact distribution annihilated by the one-forms 

Oi ~ dy — pdx, 02 ~ dp — qdx 

and that are also annihilated by the one-form 

uj = dq — F{x, y,p, q)dx. 

A contact transformation $ : — > is a local diffcomorphism that preserves the 
contact filtration on J^, meaning that 

$*6li=0 (mod6li), $*6'2 = (mod 6*1, 6*2). 

Contact transformations act on the set of differential equations by composition. 
The paper is concerned with structures that are invariant with respect to trans- 
formations of this form. An alternative characterization of third-order equations 
and contact transformations relying only on the first jet space and the contact 
structure associated to 9i will be given in Section 
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The equivalence problem for third-order ordinary differential equations has a 
long history. It was first studied under the class of point transformations by Elie 
Cartan |Car55| . In 1940, S. S. Chern |Che40| focused on the equivalence problem 
under contact transformations. A certain scalar invariant of the structure, called 
the Wiinschmann invariant, divides the family of third-order equations into two 
classes. Those with vanishing Wiinschmann invariant admit a natural conformal 
Lorentzian structure on the space of solutions. For Chern, this conformal struc- 
ture presented itself in the form of a normal 5*0(3, 2) Cartan connection on the 
solution space or, what is the same, a conformal Lorentzian metric on that space. 
A geometrical description of this conformal structure was presented much later by 
Fritelli, Kozamch, and Newman |FKN01j . who showed that the general third-order 
equation with vanishing Wiinschmann invariant could be obtained by considering 
one-parameter families of null hypcrsurfaces in a 3-dimensional space with a con- 
formal Lorentzian metric. More directly, as discussed in Section [TT] this case can 
be understood also in terms of the null geodesic spray on the bundle of null rays: 
null hypcrsurfaces then being related by means of an envelope construction. 

Chern also determined all of the contact invariants of the general third-order 
equation in which the Wiinschmann is nonzero. This was presented in the modern 
language of bundles and connections by Sato and Yoshikawa |SY98| , who in addition 
clarified the geometrical interpretation of the structure as a normal SO{3, 2) Cartan 
connection on the space J^. Nurowski and the second author later showed the 
existence of a conformal 0(3, 3) structure on a certain fiber bundle over (see 
|Nur05j ). This 0(3, 3) structure is of Fcfferman jFef76| type, in the sense of Graham 
jGra87| . if and only if the Wiinschmann invariant vanishes. Godlihski |God08j then 
proved that the associated normal 5*0(4, 4) Cartan connection included the Chern- 
Sato-Yoshikawa connection as its o(3, 2) part. 

The structure of Chern-Sato- Yoshikawa can doubtless be understood directly in 
terms of the geometry of the space of solutions of the differential equation. However, 
whereas when the Wiinschmann invariant vanishes, there is a standard geometry 
underlying the presence of certain connections — namely a conformal Lorentzian 
metric on the space of solutions — when the Wiinschmann invariant is nonzero, 
such an underlying geometry appears to be missing. The main task of the present 
article is to supply the missing geometry and to examine its precise relationship 
with these constructions. 

Invariantly associated to the structure on is a degenerate conformal Lorentzian 
metric (see Nurowski |Nur05j ) whose degeneracy is in the direction of the total de- 
rivative vector field coming from the differential equation. The second task of the 
present article is to show that this degenerate conformal structure arises naturally 
from elementary geometrical constructions on a certain curve (the polar curve) in 
the projective cotangent bundle of the space of solutions. The Wiinschmann in- 
variant itself is precisely the projective curvature of this curve. Finally, the paper 
introduces a new conformal invariant of fourth order in the metric that gives a com- 
plete geometrical characterization of degenerate metrics that arise in this manner 
from third-order differential equations. 

The first innovation of the paper is that of an incidence relation on the space 
of solutions, described in Section 13.11 An infinitesimal or linearized version of this 
idea is implicit in Wiinschmann's |Wiin05| investigations into Monge equations of 
the second degree (see also |Lie05j and [Che40| ). When the Wiinschmann invariant 
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vanishes, two solutions are incident if and only if they lie on the same null geodesic. 
When the Wiinschmann is nonzero, the incidence relation still defines a decent 
structure in a sense that is axiomatized in Section |31 roughly, the sheaf of curves 
defining the incidence relation is envelope-forming. Such a family of envelope- 
forming curves is dubbed a causal geometry, the terminology suggested by an affinity 
with structures that typically arise in the study of hyperbolic partial differential 
equations. The space of incidence curves in the solution space corresponds naturally 
to the points of the 1-jet space J^. When the incidence curves through a point are 
linearized at that point, the resulting cone in the tangent space resembles the null 
cone associated to a Lorcntzian structure. The null cone projects to a curve in the 
projective tangent space, called the indicatrix curve, borrowing terminology from 
optics |Arn97| . 

The indicatrix gives rise to a Lagrangian in a natural manner that can be writ- 
ten down in terms of the general solution of the differential equation, as described 
in Section 13.21 The Lagrangian is a function on the tangent bundle which is ho- 
mogeneous of degree two with respect to the scalar homothcty of the bundle. It 
is not fully contact-invariant, but its locus of zeros in the projective tangent bun- 
dle is invariant, and coincides with the indicatrix. Null geodesies — extremals of 
the Lagrangian along which the Lagrangian vanishes identically — are precisely the 
incidence curves. The resulting structure is a Finsler |Finl8| analog of conformal 
Lorcntzian geometry. 

The Lagrangian is in addition regular at every point of the indicatrix, and there- 
fore gives rise to a Hamiltonian on the cotangent bundle, which is described in 
Section 13.31 The zero locus of the Hamiltonian inside the projective cotangent 
space is the polar curve of the indicatrix. The total space of the indicatrix or its 
polar curve defines a 4-dimensional bundle over the space of solutions, and the 
projective Hamiltonian spray defines a projective vector field on this bundle. The 
3-dimensional quotient space under the flow of the vector field inherits a natural 
contact form from the cotangent bundle. The resulting space is contactomorphic 
to J^, the space of 1-jets in the plane, and on it the polar curves descend to a path 
geometry that deflnes a third-order differential equation in the manner described 
in Section [21 

The entire procedure sketched here is summarized in the theorem. 

Theorem 1. There is a natural local isomorphism between the set of third-order 
equations under contact equivalence and the set of isomorphism classes of causal 
geometries. 

The rest of the paper is devoted to studying the degenerate rank three conformal 
Lorcntzian metrics g on a four-dimensional space N. The fundamental invariant is 

r = g A ^v9 A ^^g A ^^g A ^^g G A^^^ ker(Fj) = S^{TN/V) 

where V is the degenerate direction. Here ker(Fj) is the space of one-forms an- 
nihilated by V, and each of the Lie derivatives ^yg lies in the symmetric square 
S^kej:{Vj). The space S^kci{Vj) is six-dimensional, and its fifth exterior power 
is naturally isomorphic to the symmetric square S'^(TN/V) of the quotient of the 
tangent bundle of by the vertical direction V. 

Theorem 2. The degenerate Lorentzian metric g on the four-manifold N arises 
(locally) from a third-order differential equation if and only if either 
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(1) T is nonzero and the classical adjoint ofT vanishes identically (equivalently, 
r has rank 1). In this case, there exists a natural conformal isometry of 
N with equipped with its invariant degenerate metric coming from a 
third-order differential equation with non-zero Wiinschmann invariant. 

(2) -2V<7 is proportional to g. In this case, the Wiinschmann invariant van- 
ishes and there exists a conformal isometry of N with equipped with its 
invariant degenerate metric coming from a third- order differential equation 
that is natural up to a gauge transformation of . 

2. Third-order differential equations 

A third-order differential equation under contact equivalence can be conveniently 
regarded as the following data. 

(1) A three-dimensional contact manifold J^. 

(2) A generic three-parameter family of (unparameterized) contact curves in 
J\ 

The contact structure on can represented as a one-form 9 which is only invari- 
antly defined up to scale. By Darboux' theorem, there exist coordinates {x, y,p) on 

such that 6 = dy — pdx. A contact curve is a curve whose tangent annihilates 
9 at each point. The family of curves is generic if at each given point of and 
tangent direction v at x annihilating 9 there exists a unique curve through x with 
tangent along v. These curves are identified with the (prolongation of) solutions 
of the differential equation. Two such structures are locally equivalent if there is a 
local diffeomorphism of to itself that preserves the contact structure and sends 
one system of curves to the other. 

Given a third-order differential equation, it is clear how to generate such a struc- 
ture by prolongation (see, for instance, [01v95| ). and the resulting structure depends 
only on the contact-equivalence class of the differential equation, by Backlund's 
theorem. Conversely, suppose we have chosen coordinates {x,y,p) on such that 
9 = dy — pdx. The distinguished class of curves is of the form 

X = x{s;a,b,c) 

y = '4'{s;a,b,c) 

p = 7r(s; a, 6, c) 

where a, 6, c are the three parameters defining a curve in the class, and s is an 
evolution parameter of the curve. There is a gauge freedom in selecting the pa- 
rameterization s of the curve, and so this freedom is eliminated this by imposing 
the condition dx/ds = 1 (that is, by effectively taking x itself to be the parame- 
ter). The contact relation takes the form p = dy/dx. Imposing this relation and 
differentiating tjj three times gives the system of equations 

y = ip{x;a,b, c) 

y' = ipx{x;a,b,c) 

y" = ipxxix;a,b,c) 

y'" = ijjxxx{x]a,b, c). 

Solving the first three equations for a,b,c in terms of x,y' ,y" and substituting 
the result into the third equation gives a third-order differential equation y'" = 
F{x,y,y',y"). 
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Associated to the jet space there is naturally a fibration tt : — > J^. This is 
the subbundlc of the projective tangent bundle PT of given as the annihilator 
of 6: =9^. Specifically, is given fiberwise by 

= {v e ¥T^J'^ I vj6 = 0} . 

It is a four-dimensional space fibered over with 5^ fibers. The space defined 
here supports the following contact-invariant structure, independently of the differ- 
ential equation. This characterization is the four-dimensional analog of structures 
studied in five dimensions by the authors in |DHS09j . 

Lemma 1. (1) There exists a natural filtration 

of the tangent bundle of . Here is the vertical distribution for the 
fibration — > J^ , is a tautological bundle of 2-planes, and is the 
annihilator of the pullback of 9. 
(2) On sections. [T{T'^),V{T^)] = [T{T^),T{T^)] ^T{T^) and[T{T^).T{T^)] = 
r(T4). 

Specifically, is the tautological 2-plane bundle whose fiber at a point (x, u) S 
Jj C PTj; J-'^ consists of all vectors v such that 7r,w is in the direction of u. In terms 
of the {x,y,p) coordinates on J^, any vector field of the form d/dx+pd/dy + qd/dp 
annihilates the contact form 9 ~ dy~pdx. Therefore this q defines a fiber coordinate 
for the fibration — >■ that allows the vector fields generating T^J^ to be 
expressed as X = d/dq, the vertical vector field for the fibration, and d/dx + 
pd/dy + qd/dp. The lemma follows by taking commutators. 

Lemma 2. Any two four-manifolds equipped with this structure are locally iso- 
morphic: a filtration C C C on the tangent bundle, such that 

[r(Ti),r(r2)] = [t{t^),t{t^)]^t{t^) and [r(r3),r(T3)] ==r(T4). 

Proof. Let Af be a manifold equipped with such a filtration and let be a nonva- 
nishing one-form annihilating T^M. Let N be the 3-manifold obtained by passing 
to the (locally defined) quotient modulo the flow of T^M. The distribution T^M 
is Lie derived along T^M, and so descends to a distribution of 2-planes on N. 
Nowhere is this distribution Frobenius integrable, and so it defines a contact struc- 
ture on N. By Darboux' theorem, N is locally contactomorphic to with its 
standard contact structure and coordinates (x,y,p). Letting g be a fiber coordi- 
nate on M N , X = d / dq generates T^M and 9 = dy — pdx. Plane subbundles 
of on which .^x maps surjectively onto are all related by a change in the 
fiber coordinate q. Indeed, as X commutes with d/dp, the latter vector field does 
not lie in T^. The bundle must contain a solution Y of .SfjcF = d/dp. One 
such solution is F = d/dx -\- pd/dy -\- qd/dp, and the ambiguity, modulo T^, in the 
solution is a transformation of the form Y i-^ Y + Xd/dp where X{X) — 0. Noting 
that A is independent of q, this ambiguity in the choice of Y can be absorbed into 
a change of coordinates q t-^ q-\- fj, where fi(x, y,p) satisfies the differential equation 
^ Ma: ^ /^y ^ Mp ~ ^- So we are free to choose the fiber coordinate q so that is 
generated by X and Y = d/dx -t- pd/dy -t- qd/dp. The coordinates {x, y,p, q) now 
defined on M establish a local diffeomorphism with ,P that sends T^M, T'^M , and 
T^M to J2, T2j2, and T'^J^, respectively. □ 
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The differential equation is specified in terms of a splitting of the first level of 
the filtration T^,P C T^J^. This is achieved by means of a vector field V of 
the form V = d/dx + pd/dy + qd/dp + F{x,y,p,q)d/dq representing the total 
derivative. While V itself is not contact-invariant, the splitting direction consisting 
of all multiples of V is. This splitting can invariantly be described in terms of 
the system of curves on that gives the differential equation. Lying over a point 
X G J^, the point u G oi the fiber is by definition a projective tangent vector 
at X. Passing through x in the direction defined by u is a distinguished curve of 
the differential equation. This curve determines a tangent direction at every point 
which therefore specifies a lift to J^. The vector V(^x.u) is the tangent direction to 
the lifted curve at the point {x, u) E J'^ . 

2.1. Conformal structure. It is possible to construct from these data a degener- 
ate conformal Lorentzian metric g on. J^. The degenerate direction for the metric 
is V, and the vector field X is null with respect to the metric. In terms of the 
coordinates {x, y,p, q) on J^, this metric is given by 

(1) g = 2[dy - p dx] [dq -\Fqdp + K dy + {\qFq - F - pK)dx] - [dp - q dx]^ 

where K — ^V{Fq) — ^F^ — ^^p- Section [SI the present paper constructs this 
metric in a manifestly contact-invariant fashion. The (locally defined) quotient 
space M = J^/V is a three-manifold, but the metric g may not be Lie derived up 
to scale along V, and so need not pass down to the quotient. The Lie derivative 
is proportional to g if and only if the Wiinschmann invariant of the original 
differential equation vanishes. The Wiinschmann invariant is given by 

(2) W = Fy + iV- lF,)K. 

Vanishing of this invariant is a necessary and sufficient condition for M to possess 
an invariant Lorentzian conformal structure. 

Going the other way, let M be a conformal Lorentzian 3-manifold. Define S to be 
the bundle over M with fiber that, at each point P, consists of all null directions 
in PTpM. The puUback metric is degenerate in the vertical direction, and so S 
supports the structure of a degenerate conformal Lorentzian 4-manifold for which 
the degenerate direction is a conformal Killing symmetry. There is a canonical 
symplectic potential "ip defined on the total space of the cotangent bundle of M. 
The form ^ is annihilated by the scaling in the fiber, and is Lie derived up to scale, 
and so descends to give a form 6 up to scale on S. The condition 9 Ad9 ^ follows 
since A dip on T*M does not vanish when pulled back to nonzero sections M — ^ 
r*M. Indeed, in local coordinates x = {x^,x'^,x^) on M with fiber coordinates 
P = iPi,P2,P3) on r*M, ip hdip = (p • dx)(dp • dx) = — ^(p x dp) • (dx x dx). Since 
dx X dx has three linearly independent components, and p x dp vanishes only on 
vectors parallel to the generator of scalings in T*M, ip /\ dip does not vanish when 
pulled back along any section of the projective cotangent bundle, and so a fortiori 
it does not vanish when pulled back to S. 

The vector field X is the given by the null geodesic spray in §. To describe this, 
fix a metric g in the conformal class on M. On r*M the geodesic Hamiltonian 
is H = T^* g^^{ip,ip), which gives rise to the Hamiltonian vector field iJ, defined 
by H_idip = dH . An integral curve )jl oi X projects to a geodesic of M, and the 
fiber component of ^jl is the covelocity of the geodesic. The image of S under the 
map S — )■ r*M is the null cone at every point of M. This is everywhere tangent 
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to the spray H, because a geodesic being initially null will always remain null. 
Furthermore, H scales quadratically in the cotangent bundle, and descends to a 
direction field on the projective cotangent bundle PT*M. This direction field is 
everywhere tangent to S, and so restricts to a direction field X on §. 

Passing to the (locally defined) quotient by X gives the space of null geodesies — 
the twistor space T. It follows from the definition of X and 9 that XaO = and 
X_id6 — 0. The second assertion follows by pulling back Xjdip — dH to the null 
cone, and using the fact that H vanishes identically there. Thus 9 descends to 
a contact structure on the twistor space. The fibers of § ^ M project to the 
distinguished curves of the twistor space. 

To prove that these distinguished curves are suitably generic, we verify that § 
admits a structure satisfying the conditions of Lemma |21 Let be the bundle 
spanned by the null geodesic spray X and be annihilator of 9. From 9 Ad9 ^ 0, 
it follows that is not Frobenius integrable, and so [r(r3), r(T3)] = r(TS). 
Moreover, since 9 is annihilated by X and Lie derived along it, [r(T^), r(r'^)] = 

r(r3). 

It remains to identify and to show that [r{T^),r{T'^)] = T{T^). Let V be 
a nonvanishing vertical vector field for § — !> M and let = span{X, V}. One 
such vector field V can be given in terms of the angular momentum operator 
L = p X d/dp. Then L is tangent to the null cone, because it annihilates H. 
On homogeneous functions of degree 0, the angular momentum factors through 
a scalar operator L(/) = V{f)p, which defines the vector field V. Because 9 
is horizontal, any such vector field, being vertical, annihilates 9. Furthermore, 
[X,V]^9 = ^x{Va9) = as well, so [X,V] e T^. It remains only to show that 
X, V, [X, y] are linearly independent. Each of these vector fields is at most first 
order in the metric, and therefore independence follows by a calculation in normal 
coordinates. In these coordinates, X = — p ■ 9/9x and so on the one hand 

[X,L]=px ^ 



X 



and on the other hand, on the null cone this acts as p[X, V\ on functions homoge- 
neous of degree zero. Finally, 

d0([X,L],L) = (p-p).9-p®p 

which vanishes nowhere on the null cone. Thus [r(r^), r(T^)] = r(r'^). 

So modulo the explicit construction of the degenerate metric and the assertion 
that specifically it is the Wiinschmann invariant that governs whether the degen- 
erate metric descends to the 3-manifold, we have proven the following theorem 
(proven in Fritelli, Kozameh, Newman |FKN01| by entirely different methods): 

Theorem 3. There is a natural local equivalence between third-order differential 
equations under contact transformations with vanishing Wiinschmann invariant and 
conformal Lorentzian 3-manifolds. 



3. Causal geometries on the space of solutions 

This section develops the natural geometric structure associated to the space of 
solutions M to a third-order equation. This structure reduces to a conventional 
Lorentzian conformal structure if and only if the Wiinschmann invariant vanishes. 
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To motivate this discussion, the previous section estabhshes the basic properties of 
the (possibly locaUy defined) double fibration 




An individual fiber of the submersion S — >■ M projects down to give a trajectory 
solving the differential equation in T: this is a re-expression of the notion that M is 
a space of solutions of the differential equation. The fibers for the other submersion 
§ — !> T also project down to the space of solutions M, although their precise meaning 
has heretofore not been identified in general. 

When the Wiinschmann invariant vanishes, M carries a natural conformal Lorentzian 
metric by Theorem [31 and S is canonically identified with the null cone bundle asso- 
ciated to this metric. The space T is then the twistor space: the quotient of § by the 
null geodesic spray. The fibration S — >■ M can be understood as the subbundle of 
the projective tangent bundle PTM of null directions. Under this correspondence, 
a null cone with vertex at P G M corresponds to a one-parameter family of null 
geodesies, which in turn is identified with the trajectory defining the solution P. 

This structure can be axiomatized in a manner that allows construction of the 
spaces S and T, along with their natural contact structure. The one parameter 
families of null geodesies starting at each point P give rise to a cone in the tangent 
space TpM with vertex at the origin. Equivalently, such a cone is the affine cone 
over some curve in the projective space PTpM. Here locality considerations may 
mean that the cone may fail to close up completely, or the associated curve may 
have one or more singular points. Henceforth, we shall work only near regular 
points of the curveQ 

3.1. Incidence relation and the indicatrix. Two solutions yi and 1/2 to the 

differential equation y'" = F{x, y, y', y") are said to be incident if, at some point a;, 
one has 

yi{x) = y2{x), y[{x) = y'2{x). 
That is to say, two curves are incident if and only if the associated solution curves 
in the intersect. In the latter interpretation, the incidence relation is manifestly 
contact-invariant . 

The set of all solutions incident with a given solution P G M cuts out a surface 
Np with a conical type singularity at the vertex P. The generators of this cone 
can be described as follows. If Q is a solution incident with P, then the curve 
NpQ consisting of solutions R incident with P at the same point of as Q. (Our 
localization assumption implies that two distinct solutions are incident at most at 
a single point of J^.) The surface Np is thus ruled by the pencil of curves NpQ as 
Q varies over solutions incident with P. 

The incidence relation described here is called a causal geometry. The following 
axioms are assumed to hold in a sufficiently small neighborhood of each point 
P G M: 



Wc can localize near a point of T and a point of M simultaneously. Thus all solutions can 
be assumed to be fully regular, but cannot necessarily be continued for all time. The geometrical 
implication is that we may isolate a smooth part of the null cone at each point of M, but the 
"cone" need not then close up. 
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Axiom 1. Np \ {P} is a smooth embedded hypersm'face in M 

Axiom 2. The singular variety Np is ruled by a pencil of smooth curves NpQ from 
P to Q, as Q varies over Np \ {P}. 

Axiom 3. For every Q G Np\{P}, the curve NpQ on Np coincides with the curve 
Nqp on Nq: Npq ^ Nqp. 

Axiom 4. Let Cp be the set of tangents at P to the curves Npq as Q varies 
over Np \ {P}- Then Cp is a regular (smooth) curve in the projective 
tangent space PTpM, and does not make second order contact with any 
of its tangent lines. (The curve Cp is called the indicatrix.) 

Axiom 5. Let Q,R ^ Np \ {P} be given distinct points which are mutually inci- 
dent. Then the tangent plane to Np at R is the same as the tangent 
plane to Nq at R: TpNp — T^Nq. (Envelope condition) 

Theorem [5] establishes that these axioms hold in the case of the causal geometry on 
the space of solutions of a third-order differential equation. The first four axioms 
are fairly natural assumptions that make precise the notion that the Np should 
be a cone based at P. The envelope condition of axiom 5, in terms of the bundle 

— > M, guarantees that the one-form 6, whose annihilator is the 3-plane bundle 
that lifts the tangent planes through the vertices of the cones Np, is Lie derived 
along the fibers of the fibration ^ . 

Causal geometries also appear naturally in connection with the following varia- 
tional problem. Let L : TM \ {0} — K. be a homogeneous function of degree two 
(a Lagrangian), possibly only defined on an open conical subbundle of TM \ {0}. 
Suppose that for each P G M, the curve Cp = {w G PTpM | L{P, v) = 0} is smooth. 
For a curve 7 in M, consider the energy functional 



Eb] = l I L{^{t),i{t))dt. 

A null geodesic is a curve 7 that is an extremals for E along which L{'-f{t),j' (t)) = 0. 
The family of null geodesies for a Lagrangian defines a causal geometry if and only 
if the Lagrangian is regular at every point of the cone over Cp in TpM for each P, 
in the sense that its vertical Hessian is nondegenerate in directions tangent to Cp. 
In coordinates a;* for M and fiber coordinates for TM, this is the condition that 
the 3x3 matrix a-ia - , be nonsingular in directions tangent to the null cone. This is 
also a sufficient condition for L to be conserved along an extremal, and thus the null 
geodesies are precisely those extremals of the energy for which ^(7(0), 7'(0)) = 0. 

The Euler-Lagrange system is a second order ordinary differential equation for 
the curve "f(t). By smooth dependence on initial conditions, for initial conditions 
lying on the hypersurfacc L{P,j'{0)) — sufficiently near the origin of the tan- 
gent space TpM, the null geodesies foliate a smooth hypersurfacc in M. giving 
axioms 1-2. Axiom 3 follows since the null geodesies are critical points of the en- 
ergy under compactly supported variations, and so in particular are characterized 
independently of the direction of their parameterization. Axiom 4 follows from the 
regularity of the Lagrangian. Finally, for axiom 5, it is sufficient to show that, for 
any point P, and any variation 7^ of null geodesies through P, 



(3) 



= 0. 

s=0 
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This then establishes that the nuh cone at P is tangent to the indicatrix curve at 
every point, which is equivalent to axiom 5. By the Euler-Lagrange equations, 

= 1l(7, 7,) = — ^ + — ^ 
ds ' dx"^ ds dx^ ds 

d dL \ dY dL dj' 
dt dx^ J ds dx^ ds 
d f dL d-y' 



dt ydx"^ ds 

So the left-hand side of @ is constant along the curve 7(t), as required. 

Conversely, a Lagrangian can be associated to any causal geometry, up to a 
certain ambiguity. Let L : TM ^ R be a function homogeneous of degree two that 
vanishes at every point of the cone over the indicatrix curve Cp for every point 
P e M, and suppose further that the vertical Hessian of L is nondegenerate. ft 
is always possible to select such a function, at least locally, and the ambiguity in 
selecting such a function is of the form 

L{x, x) fl{x,x)L{x,x) 

where fi is a nonvanishing function. A causal geometry is already associated to such 
a Lagrangian, by the preceding argument, and so it is sufficient to show that the 
causal geometry induced by the Lagrangian agrees with the one already given. To 
show that the curves Np coincide with the null geodesies of L, it is enough to show 
that they are critical points for the energy with respect to variations 7s tangent to 
the null cone, meaning that holds. Indeed, the Euler-Lagrange system for the 
constrained problem is 

dL d dL _^dL 

dx^ dt dx"^ 
where A = A(t) is a free function of t. The freedom in A can be absorbed by a 
change in parameterization of the curve '){t). 

Now consider the curve 7 = Npg, with some parameterization. Let 7^ be a 
variation of 7 obtained by varying the endpoint Q along Np (that is, 7^ is a vari- 
ation through null geodesies). In the statement of Axiom 5, taking the point R to 
approach Q, a necessary condition for 7^ to be on the null cone is that 1^ holds. 
Moreover, also by Axiom 5, 7^ satisfies L{'ys,js) = 0. Therefore we have on the 
one hand 

Q d . ^ dL c?7* ^ dL 7* 

ds ' dx^ ds dx'- ds 

And on the other hand also 

d ( dL df\ f d dL\ d-y^ dL df 



dt \ dx'' ds J \dt dx' J ds ' dx'' ds 
Equating the two gives 

dL d dL\ dj 



dx' dt dx' J ds 
along NpQ. Thus Npq is a null geodesic for the Lagrangian L. 
These considerations are summarized in the following theorem: 

Theorem 4. Let L : TM R be a homogeneous function of degree two, such that 
for each P G M, the curve Cp = {L = 0} C PTpM is smooth. Suppose also that the 
vertical Hessian of L is nondegenerate at each point of Cp. Then the null geodesic 
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curves of L form a causal geometry. Every causal geometry arises (locally) in this 
manner. 

3.2. The causal geometry associated to a differential equation. At each 
point P of M, denote by 

y = f[x-P) 

the solution to the differential equation represented by P. We break contact invari- 
ance in specifying the differential equation, but shall ultimately be concerned only 
with the contact-invariant information that can be extracted from /. Locally, this 
solution depends smoothly on x and P. The underlying assumption under which 
M is a differentiable manifold smoothly parameterizing a space of solutions is that 
in local coordinates P = (pi,p2,P3) on M, the Wronskian determinant 



(4) 









dpi 




dp3 


dfjn 






dpi 


dp2 


dps 








dpi 


dp2 


dps 



7^0. 



To obtain a more concrete description of the cone Np through a particular 
solution P S M and the associated curves NpQ that rule the cone defined in the 
previous section, suppose that j(t) is a parameterization of the curve NpQ such 
that 7(0) = P. Since all points along 7(t) are incident with P at the same point 
(x,?/, y'), the following two equations must hold along 7: 



(5) 



/(x; 7(0) = /(a;; 7(0)) 
/,(x;7(t)) = /x(x;7(0)). 



Under generic conditions, the second equation can be used to solve for x in terms of 
7(0) and 7(t), and then substituted into the first equation which may then be solved 
for the admissible values j{t). Although ostensibly this system of equations is not 
invariant under contact transformations, by Backlund's theorem a contact trans- 
formation will preserve the space of solutions j{t). This is geometrically evident 
because the causal geometry itself is contact-invariant. 

Solutions of these equations arise as first integrals of the differentiated forms 



(6) 



dt 
d_ 
dt 



fix-nit)) ^0 
/,(x;7(i)) = 



with initial conditions 7(0),7'(0). The initial conditions are not completely arbi- 
trary. Rather if 7(0) is fixed, then 7'(0) is constrained by the requirement that 



(7) 



d 
dt 
d 
dt 



fix-.lit)) 
Ux-nit)) 



= 



t=o 



= 



t=o 



The second equation of ([7]) can be solved to obtain the parameter x in terms of 
the initial conditions 7(0), 7'(0), by (0]). The following Lagrangian is homogeneous 
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of degree two in 7'(0): 

(8) i(7(0),7'(0)) - -^/(x(7(0),7'(0));7(t))|/..(:^(7(0),7'(0));7(t)) 

at at 

The second factor ensures that the Lagrangian is regular (Lemma [3]) . For fixed 
7(0), the values of the tangent 7'(0) satisfying equation ([7]) cut out an affine cone 
over a curve in the projective tangent space at 7(0). This cone is the "null cone" 
for the Lagrangian L, and it coincides with the tangent cone to N^{o) at the vertex 
7(0). As in Axiom 4, for P e M, denote by Cp C PTpM the curve cut out in the 
projective tangent space by the equation L(P, v) = 0. 

Although the Lagrangian L is not itself contact-invariant, by the argument al- 
ready given its locus of zeros is contact-invariant. Under contact transformations, 
the Lagrangian is determined up to rescaling by a nonvanishing function of 7(0) 
and 7'(0), 

L(7(0),7'(0)) ^ 17(7(0), 7'(0))L(7(0), 7' (0)). 

The curves Npg that generate the causal cone Np are extremals for the energy 
functional 

m^l J'mt),i'it))dt. 

Theorem 5. A third-order differential equation determines a causal geometry sat- 
isfying Axioms 1-5. 

Proof. Axiom 1 reflects our running localization assumption that the third-order 
differential equation under consideration is regular. In local coordinates on M, the 
system of equations 

(9) f{x;qi,q2,q3) ^ f{x;pi,p2,P3), fx{x;qi,q2,q3) = fx{x;pi,P2,P3) 

for unknowns x,Q = ((Ji, 92j 93) and fixed P = {pi,P2iP3) has Jacobian matrix 

fxix;qi,q2,q3) - fxix]Pi,P2,P3) ^ ^ ^ 

fxx{x;quq2,q3) ~ fxx{x;pi,P2,P3) ^ ^ ^ 

The matrix always has rank two by (jlj. Moreover, at a solution Q ^ P oi the orig- 
inal system, the lower left-hand corner cannot be zero, by uniqueness of solutions. 
Thus the first column, along with one of the remaining three columns must yield 
an invertible 2x2 submatrix. The implicit function theorem then implies that the 
solution is a smooth submanifold away from the vertex Q ~ P. 

For axiom 2, Npp is the set of solutions Q = (gi, 92, 93) of ([9]) for a given value 
of x,P = (pi,P2,P3)- Because ([9|) has rank two in the qt variables, the space of 
solutions is a smoothly embedded curve. These curves clearly cover Np. Further- 
more, two distinct curves meet only at the vertex P, by uniqueness of solutions of 
differential equations. 

Axiom 3 is obvious. Axiom 4 is equivalent to the assertion that L is a regular 
Lagrangian, which is proven in Lemma |3] of the next section. 

Finally, suppose without loss of generality that in the statement of axiom 5, R is 
a point between P and Q on NpQ . Along the curve NpQ the value of the parameter 
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X is fixed, and Npq itself consists of all points R such that 

Ux-R) = Ux-P). 

The tangent plane to Np at the point R is the annihilator of di{J{x,R)^ This is 
the same tangent plane as that obtained by interchanging the roles of P and Q. □ 

3.3. Associated Hamiltonian. As in the previous section, for P G M, let y = 

f{x;P) denote the solution of the differential equation corresponding to P. Once 
again, contact-invariance is broken, but ultimately we will only be concerned with 
contact-invariant information contained in the solution. Let dpf denote the exterior 
derivative of / regarding x as constant. In local coordinates (pi,p2,P3) on M, 

dpj(x;pi,p2,P3) = i^dpi + -7—dp2 + j^dps. 

dpi dp2 op3 

Then x i— ?■ dpf{x;P) defines a curve in the cotangent space TpM.. Denote the 
associated projective curve by Cp C Pr*M. 

This curve is linked to the curve Cp cut out by the Lagrangian via the following 
construction. Let ^ be a three-dimensional vector space and C a smooth curve 
in the projective plane FV. The dual projective plane FV* is naturally identified 
with the space of lines in FV. The dual curve C* is the curve in FV* defined by 
locus of lines tangent to C. Suppose that C is a nondegenerate curve in FV, with 
parameterization t i— > j{t). At a point j(t) of C, the corresponding point of the 
dual curve is obtained by solving for 7*(t) G FV* the equations 

(7*(t),7W)=0 

(7*W,7'W)=0. 

Properly speaking, to make sense of the second equation, it is necessary to choose 
a lift of 7 to a curve in V. Modulo the first equation, the second equation does 
not depend on the choice of lift, and so there is no ambiguity in speaking of the 
solution of the system of equations. 

Proposition 1. Cp C Pr|,M and Cp C PTpM are mutually dual. 

Proof. The curve Cp is characterized as the image of the map 

X H- dpf{x; P). 

The dual curve to Cp is defined by the equations 

{l*{x),dpf{x-P))=Q 
^ ' {j*{x),dpUx;P))^0. 

But these two equations are identical with the equations that characterize Cp. □ 

Alternatively, choose local coordinates at P and linearize the differential equation 
at the solution defined by P. Then 

f{x;pi,p2,P3) = (f>i{x)pi + (I)2{x)p2 + 03(x)p3 

and so 

dpf = (j)i{x)dpi + (I)2{x)dp2 + (t)3.{x)dp'i. 



^Here and elsewhere, the notation dfif{x; R) is the exterior derivative of / with respect to the 
variable R only. Equivalently, it is the exterior derivative of / modulo the relation dx = 0. 
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The incidence relation between f{x]pi,p2,P3) and a nearby solution f(x;pi + 
dpiiP2 + dp2,P3 + dps) is then precisely 

(j)i{x)dpi + cj)2{x)dp2 + (f)3{x)dp3 = 

(j)[{x)dpi + (j)'2{x)dp2 + (t)'y,{x)dp3 = 

but these are the same equations that characterize the dual curve of Cp. 

Lemma 3. Let L : TM — > M 6e a Lagrangian (jlip . Then L is regular in a 
neighborhood of each point of Cp. 

Indeed, it is sufficient to show that the Hessian matrix L / dqidqj is nonsingular 
at each point of Cp. The Lagrangian is defined by 

(11) L{q,q) = {qi^{x{q,q),q)^ {i3-^{^{l^i)^l) 
The function x{q, q) is defined by 

(12) q,^{x{qA),q) = Q- 

In particular, by implicit differentiation, 

dx _ df^./dq.i 



^"^^^ dqt qkdfxx/dqk 

At a point of Cp, in addition the following holds: 

df 

(14) q,^{x{q,q).q)=Q- 

The Hessian at a point of Cp is computed by differentiating imposing 
((T^ and (fT?)l along the way. and then finally substituting p^ . In detail, denote 
= 9//%, /j;, dfx/dqi, etc., and = dx/dq^. Then 

= (tfk(ffxxl 

Li = fiq^fxxi + q'^fkfxxi + ^ fxkXi^i fxxi + q'' fkq^ fxxxiXi 
Lij fxi'^jq fxxi ^~ fifxxj ~i~ fiq fxxxH'Xj fjfxxi~^ 

fxj^iq fxxi ^ fxxk-^i-^j (ifxxi + fjq'^fxxxe.Xi (mod (HH, dH])) 

_j- |. Iff If f ^ f f ( ^ fxxxi \ r r f q fxxxl 

— Jxijxj I Jijxxj 'y J xxi J j Jijxj \ .1, » I Jxijj I . 7, ^ 

V q fxxk ) \ q^fxxk 

by Thus the Hessian matrix of L has the form 

Hess L = -AA^ + BC^ + CB^ 
where the column vectors A^B^C are defined by 

Ai — f xi 



■I f 

i ~ J xxi p jx 

Q Jxxk 



By the hypothesis (Hj), A^B^C arc linearly independent, and so 

detHessL - -(det[yl B C]f ^ 0, 
which establishes the lemma. 
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The Legendre transformation associated to the Lagrangian L is a function I£ : 
TM T*M covering the projection onto M. Over a point P e M, L : TpM R, 
and ^ = L, the vertical exterior derivative of L. The Hamiltonian associated 
to the degree 2 homogeneous Lagrangian L is defined hy H = L o ££~^ . Here 
J^"^ is the inverse, possibly only locally defined near points of Cp, of the function 
^ : TpM T^M. 

By Lemma [21 H is well-defined in a neighborhood of the preimage of Cp under 
The Hamiltonian, where it is defined, vanishes precisely on the dual curve Cp. 
Indeed, in local coordinates, 

dL 

— -qT^ ~ fiQ fxxl + q fkfxxi + q fxkXiq ' fxxl + q fki fxxxt 

— fiq fxxl 

when evaluated at any point of Cp. Thus is proportional to dpf at each point 
of Cp. Inverting, we conclude that L o .^"^{dpf) = 0, so H vanishes along Cp. 

The following alternative construction of the Hamiltonian also applies, by lin- 
earizing the problem at P. In local coordinates at P, (j)i{x) = i = 1, 2, 3 define 
independent solutions of the linearized ordinary differential equation. In terms of 
these three solutions, the linearized equation itself can be recovered by solving the 
3x3 system for the unknown coefficients ht: 



(15) Mx)ho{x) + <l>-{x)hi{x) + <l>Ux)h2{x) 

Cramer's rule gives 



(bTix), z = l,2,3. 



ho{x) 











0'l 












<P'2 






















^1 
























03 


<P'3 





hi{x) 



"1 



0i 

02 



h2{x) 



The Lagrangian of the original equation localized at the point P is equal to the 
Lagrangian of the linearized equation. It is given first by solving 

(t>'i{x)qi + (j)'2{x)q2 + 03(a;)g3 = 

for function of (71,92,93- In that case, 

L{q) = (0i(a;(g))gi + 4>2{x{q))q2 + h{x{q))q^) {(l>'i{x{q))qi + (I)2{x{q))q2 + 4>'i{x{q))q:i) . 

The associated Hamiltonian is obtained by the same construction, but applied 
to solutions 01,02, 03 of the adjoint equation to (fTC)) : 

y{x)ho{x) - {y{x)h^{x))' + {y{x)h2{x))" = -y"\x). 

The following theorem is due to Wilczynski [Wil06| ; cf . also Olver [01v95) : 
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Theorem 6. The projective curves x i—)- 4'i{x) d/dqi + 4)2{x) d/dq2 + (pzi^) d/dq^ 
in PTpM and x H> 0i(a;) dqi + 02 (a;) dq2 + 03 (a;) dq^ are mutually dual. 

4. HAMILTONIAN SPRAY 

As the point P varies, the dual curve Cp cut out by the Hamiltonian defines a 
subfibration of the projective cotangent bundle Pr*M. The four-manifold defined 
by the total space of this fibration is denoted here by S, and the cone over S in 
the cotangent bundle r*M is denoted by S. The canonical one-form 6 on T*M is 
homogeneous of degree one, and pulls back to a natural one-form defined up to 
scale on §. 

For a fixed choice of Hamiltonian homogeneous of degree two. define the Hamil- 
tonian spray on r*M as the unique vector field X such that 

X_ide ^ dH. 

The Hamiltonian spray is invariant up to scale under rescalings of H . Moreover, it 
is tangent to the variety § cut out hy H = Q since X_idH = 0. The vector field X 
is homogeneous of degree one: if fit ■ T*M T*M denotes the dilation mapping 
in the fibers, then Xta = t{fit)*Xa. 

The Hamiltonian spray will now be used to define a filtration on TS in a manner 
analogous to the proof of Theorem [31 Let y be a nonvanishing vector field that 
is vertical for the fibration § M. Let C TE> be the subbundle spanned 
by X, the subbundle spanned by X, V. Let T'^ C TS be the annihilator of 
6. Since X and V both annihilate 9, C T'^. Moreover, since 9 is annihilated 
by X and is Lie derived along it, [T{T^),r{T'^)] = T{T^). Because 9 A d9 does 
not vanish when pulled back along on any section of the cotangent bundle, it also 
does not vanish on § and so is not Frobenius integrable at any point, and thus 

[r(T3),r(r3)] = r(T§) 

It remains only to show that [r{T^),T{T^)] = r(r^). It is sufficient to prove 
that X, V, [X, V] are linearly independent. As in the proof of Theorem [31 it is 
convenient to work with a particular choice of vector field V. Define the tensor 
h e Sym^ T{T*M) to be the vertical Hessian of H. In coordinates, 

and let hij be the inverse of h'^ . These two tensors can be used to raise and lower 
indices. Let e be the associated volume tensor in the fiber. In coordinates 

e = e'^'^dpi (E) dpj ® dpk- 

The angular momentum 

Q 

V = e^^^pjliMTT- 
ope 

kills H, since OH /dpi = p^, and so is tangent to the null cone bundle S. Fur- 
thermore, on homogeneous functions of degree zero along the null cone, factors 
through a scalar operator = p^V, because p* and d/dpi are an orthogonal basis 
for the orthogonal complement oipi. This defines a vertical vector V . To show that 
X, V, [X, V] are linearly independent, it is enough to show that d9{[X, U],U) ^ 0. 
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The Lie bracket is given by 
for some tensor T. So 

which does not vanish when restricted to the null cone. 

Thus by Lemma [5J S equipped with its geodesic spray and vertical vector field 
is locally isomorphic to the space J^. Because of the preferred direction V , the 
inclusion of bundles C splits, and thus gives rise to a third-order differential 
equation the distinguished curves of which are the fibers of S — M. The entire pro- 
cedure is reversible, by the construction of the preceding section, which establishes 
Theorem [TJ 



5. Recovering the degenerate metric 

The degenerate metric on S is defined as follows. A point of § consists of a point 
P e M and v £ Cp. Now, through v 6 Cp, there is a uniquely defined osculating 
conic to Cp at v. This osculating conic in turn defines a unique conformal metric 
hp^y : T*M X r*M — !• M. A degenerate conformal Lorentzian metric is defined by 
pullback on the subspace of the cotangent bundle of § that annihilates the vertical 
direction: 

gp^y{a,l3) = /ip, t,(7r, a, 7r,/3). 

A proper degenerate conformal Lorentzian metric is obtained by dualizing In 
order to derive the formula for the metric ([T]), it is necessary to obtain explicit 
formulas for the osculating conic of a projective curve. The overall program is 
inspired by the work of Wilczynski |Wil06j . 

Suppose that (^(t) = {(j>i{t),(j)2{t),4>3{t)) parametrically specifics the homoge- 
neous coordinates of a projective curve, with dot {(f), cf)' , (f)") ^ 0.. The projective 
curve is a conic provided that there exists a 3 x 3 symmetric non-singular matrix 
A such that 

(f)^A(f) = 0. 

Supposing that is given, the task is to determine a matrix A such that at a 
given point t — to this holds to as many orders in the expansion in powers of 
t — to as possible. Since A is regarded projectively, it has 5 independent numerical 
components. These are obtained by solving the system of 5 equations linear in the 



If y is a vector space and W d V , and B is a nondegenerate bilinear form on W , then B 
gives rise to a linear isomorphism : W — ^ W . The dual (degenerate) form on V' is given by 

the mapping Tg : V — > V' /W^ W — - — > W — V V where the first is the quotient map, the 
second is the natural isomorphism, the third is the inverse of Tb, and the last is the inclusion 
map. 
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entries of A: 



(0^A0)(<o) 


^ 


(0^A0)'(<o) 







= 


^^A0)"'(to) 







= 



Once such a matrix is found, the obstruction to continuing to the fifth order is the 
derivative (0-^A0)(^)(<o), and is the projective curvature associated to the curve. 

By the first two equations, A(j){to) is proportional to the cross product 4'{to) x 
(j)'(to), and since A is taken projectively, we can fix a scale by taking 

A6(t ) = '^(^o) X '^'(^o) 

dct(0(to)>'(to),0"(to)) 

or, equivalently, (l){to)'^ A(p" (to) = f. The third and fourth equations then give, 
respectively 



0'^(<o)^</''(io) = -1 
0'^(to)A0"(to) 



'T't ^A^"(t \ ^ l dct(0,0',0"Ofa) 
3det(</),0',</)")(to)' 
The final equation now gives 

■ij."T(. \ AA."n- ^ _L AA^'Tti. \ AA"'(+ \ det(0,0',0W)(to) 
Now the coordinates of the curve (^{t) satisfy a third-order differential equation 

= ho{t)<j,{t) + hi{t)^'{t) + h2{m"it) 

where hi are given explicitly in terms of determinants of cf) and its first three deriva- 
tives as in (fT5|). The above equations reduce to 



(/.^(to)A(/.(to) 


= 


cl,^{to)A4>'{to) 


= 


^'^{to)A<j>"{to) 


= 1 


{to) Ac/,' (to) 


= -1 


0'^(to)A0"(to) 




30"^(to)^<^"(to) +40'^(io)V'(to) 





The last equation simplifies by substituting (|15p for (/>'" and then using the remain- 
ing equations to give 

3(P"^ito)A^"{to) = hlito) - /i^(to) + -Shiito). 

The obstruction ((/)^A(/))(5)(to) can now be calculated by expanding any terms in- 
volving (/)"', (/)(^', in terms of lower order and then using the above equations. 
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We find that 
(16) 

(0^A</))(5)(io) - 12/io(to) + 4/ii(to) + Ih^itof - 6h[{to) - 4/i2(<o)/i2(io) + 2h'^{to), 

which is precisely the Wiinschmann invariant for the equation (|15p . 

The matrix A obtained from this procedure is also of interest, because it gives the 
conformal Lorentzian structure. In the basis ((/>, 0', the symmetric 2-tensor A 
is given by 





"0 





1 


A = 





-1 


-lh2{to) 




1 




hl(to)-3h'2ito)+9hi{to) 
9 



When, as in section [331 the curve £ TpM is the hnearization of the sohition 
f{x;P) to the differential equation at a point P G M, then (fTSj) is the linearization 
of the differential equation at P: 

ho{x) = Fy{x, f{x; P), f^{x; P), f^x{x; P)), hi{x) = Fp{x, f{x; P), fx{x; P), f^^{x; P)), 

h2{x) = F,{x, fix; P),Ux; P), P)) 

The projective curvature obtained from (jl6p agrees with the Wiinschmann invariant 
([2]) of the original equation. Substituting in for the components of the 2-tensor A 
gives the metric 

\dp J dq dy 3 dp dq 9 \9q/ 

The dual degenerate conformal metric, defined on the full tangent bundle of §, 
agrees with ([1]). 



6. Inverse problems 

By Theorem 1, every causal geometry gives rise to a third-order differential 
equation. A more subtle inverse problem is, given a rank three degenerate conformal 
Lorentzian metric on the tangent bundle of a four manifold, when is there a causal 
geometry from which it arises? More precisely, given a one-parameter family of 
plane conies, defined by 3 x 3 nondegenerate symmetric 2-tensors A{t), when is 
there a plane curve (j){t) such that, as t ^ to, 

Interchanging t and to, a necessary and sufficient condition is that 

0(<orA(t)0(to) = o(t-to)'. 



That this basis is "nonholonomic" (i.e., nonconstant) is significant for the inverse construction, 
discussed presently. 
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So at each to, the point 0(io) must satisfy five equations 

0(to)^A(to)<^(io) = 

(PitofA'itomto) = 

cl,{tofA"{to)cf>{to) = 

(bitofA'"ito)Hto) = 

0(torAW(t„)<^(fo)=O 

In this system, the matrices A{to), A' (to), A" (to), A'" (to), A'^'^'> (to) should be re- 
garded as given, and the 3-vector (j){to) as unknown homogeneous coordinates. The 
system is clearly overdetermined: the two (projective) degrees of freedom in (/>(to) 
must satisfy five equations. Geometrically the point 0(to) niust simultaneously lie 
on five plane conies, but the intersection of more than two plane conies is generically 
empty. 

The overdetermined system gives rise to a consistency condition on A and its 
first four derivatives, which we now describe. In the generic case, we can solve the 
linear equations for a 3 x 3 symmetric matrix X 

tiA{to)X = 
tTA'{to)X = 

(17) tTA'\to)X = 

tTA"'{tQ)X = 
tTA^'^\to)X = 0. 

This can be solved uniquely for X, up to scaling, provided the system has rank five. 
The consistency condition is then that the solution X has rank one and so splits as 
an outer product 

X = cj){to)Htof- 
This happens if and only if every 2x2 minor of X vanishes 

6.1. Intermediate cases. When (fTTj) has rank one, the Wiinschmann invariant 
vanishes. When it has full rank, then it gives rise to a causal curve provided the 
2x2 minors of the solution X all vanish. A calculation done in Mathematica 
shows that, for structures coming from third-order equations, these are the only 
two possibilities: either the system has full rank (and thus nonzero Wiinschmann) 
or it has rank one (and zero Wiinschmann). 

However, a priori such a system, coming from an arbitrary degenerate conformal 
Lorentzian structure in four dimensions, can have any rank between 1 and 5. It is 
interesting to understand why these intermediate cases do not lead to causal curves. 

Rank 2. Suppose that (|17p has rank two (in an interval around Iq). Then 
^"(0 = /(*)^(^) + foi' some functions / and g. The initial matrices 

A{U)), A' (to) can be brought simultaneously into diagonal form by a transformation 
of the form 

A{tQ) ^ M^A{ta)M, A'{to) ^ M^A'{to)M. 

^Thc minors are not independent, however. The variety in P(Sym^(R^)) on which the 2x2 
minors of a symmetric 3x3 matrix vanish is the well-known Veronese surface, which is not a 
complete intersection. Locally it is the zero locus of any three minors coming from distinct rows 
and columns. 
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, x = 








X3'i_ 



Relative to this fixed initial basis, A{t) and A'{t) remain diagonal throughout the 
interval of existence. It is convenient to put 

A{t) 

The two equations of (fT7)l reduce to 

ix A{t)X = A{t) ■ X = Q 

ii A'{t)X = A'^^{t)Xii+ A'^S)X22 + A'^S)X33 

Solving: 

X = A{t) X A' it) 
up to an overall scale. Now if (j){t) is a curve solving 

(j){tf A{t)4>{t) = (t){tf A' {t)4>{t) = 

then the entries of must square to the entries oi X, so that 



±y A X A' ■ ei 



±\AxA'-e2 



±y A X A' • £3 



Differentiating gives 



Ax A' 



and so 0' is proportional to 0. Thus the range of 4>{t) is a projective point (in the 
complex sense) □ Therefore there are one or no solutions, depending on whether 
the square roots are all real. 

Rank 3 and 4- We argue indirectly that, if (fTT]) has rank 3 or 4 throughout an 
interval and 4>{t) is a solution in that interval, then 0(t) parameterizes either 
a projective point (as in the rank 2 case) or a line (which is degenerate from our 
point of view). We were unable to devise a direct argument analogous to the 
rank 2 case. In general, because (f>(t) is a three-vector, there must be a non-trivial 
linear relation between (j}{t) and its first three derivatives. This must either give a 
proper third-order equation in an interval, or else there is a linear relation between 
(pjfj)' ,(/)" . In the latter case, (j){t) does indeed parameterize a line (if it satisfies a 
second order equation) or a point (if the equation is first order). In the former 
case, the coefficients of A{t) can be given in terms of (j){t) , (/>' (t) , (j>" (t) and the 
coefficients of the third-order equation hi{t), h2{t), ft.3(t), as in the previous section. 
But, as already indicated, this implies that the system (fT7|) has rank either 1 or 5, 
a contradiction. 



Indeed, if 4>{0) and 4>'{0) are linearly independent, then <^(0) X <7!>'(0) ■ </<(t) satisfies a second 
order ode with both initial conditions zero, so (/'(O) X 4>'(0) ■ (j>{t) = which implies that <f>(i) is 
constrained to a projective line. If </'(0) and ^'(O) arc linearly dependent, then smoothness of 
dependence on initial conditions gives the result. 



22 



JONATHAN HOLLAND AND GEORGE SPARLING 



References 



[Arn97] V. I. Arnol'd, Mathematical methods of classical mechanics, volume 60 of Graduate 
Texts in Mathematics, Springer- Verlag, New York, 1997, Translated from the 1974 
Russian original by K. Vogtmann and A. Weinstein, Corrected reprint of the second 
(1989) edition. 

[BDE08] R. L. Bryant, M. Dunajski and M. Eastwood, Metrisability of two-dimensional pro- 
jective structures, (2008), arXiv/0801.0300. 

[Carlo] E. Cartan, Les systemes de Pfaff, a cinq variables et les equations aux dcrivces par- 
tielles du second ordre, Ann. Sci. Ecole Norm. Sup. (3) 27, 109-192 (1910). 

[Car24] E. Cartan, Sur les varietes a connexion projective, Bull. Soc. Math. Franco 52, 
205-241 (1924). 

[Car32] E. Cartan, Sur la geometric pseudo-conforme des hypcrsurfaces de I'espace dc deux 

variables complexes II, Ann. Scuola Norm. Sup. Pisa CI. Sci. (2) 1(4), 333-354 (1932). 
[Car33] E. Cartan, Sur la geometric pseudo-conforme des hypersurfaces de I'espace de deux 

variables complexes, Ann. Mat. Pura Appl. 11(1), 17-90 (1933). 
[Car55] E. Cartan, La geometria de las ecuaciones diferencials de tercer ordcn, in Oeuvres 

completes, Partie III, volume 2, Gauthier-Villars, Paris, 1955. 
[CG06] A. Cap and R. Cover, A holonomy characterisation of Fefferman spaces, ESI Preprint 

1875 (2006). 

[Che40] S.-S. Chern, The geometry of the differential equation y'" = F{x, y, y' , y"), Sci. Rep. 

Nat. Tsing Hua Univ. (A) 4, 97-111 (1940). 
[CM74] S.-S. Chern and J. K. Moser, ReaJ hypersurfaces in complex manifolds. Acta Math. 

133, 219-271 (1974). 

[CS05] M. Crampin and D. J. Saunders, Cartan's concept of duality for second-order ordinary 

differential equations, J. Geom. Phys. 54(2), 146-172 (2005). 
[CS07] A. Cap and K. Sagerschnig, On Nurowski's Conformal Structure Associated to a 

Generic Rank Two Distribution in Dimension Five, ESI Preprint 1963 (2007). 
[DHS09] B. Doubrov, J. Holland and G. SparUng, Space-time and G2, (2009), arXiv/0901.0543. 
[Fef76] C. Fefferman, Mongc-Ampcre equations, the Bergman kernel, and geometry of pseu- 

doconvex domains, Ann. Math. 103(20), 395-416 (1976). 
[Finl8] P. Finsler, Uber Kurvcn undFF'achen in allgemeinen Raumen, PhD thesis, Gottingen, 

1918. 

[FKNOl] S. Frittelli, C. Kozameh and E. T. Newman, DiffcrcntiaJ geometry from differential 

equations. Comm. Math. Phys. 223(2), 383-408 (2001). 
[FKN02] S. Frittelli, N. Kamran and E. T. Newman, Differential equations and conformal 

geometry, J. Geom. Phys. 43(2-3), 133-145 (2002). 
[FKNN02] S. Frittelli, C. Kozameh, E. T. Newman and P. Nurowski, Cartan normal conformal 

connections from differential equations, Classical Quantum Gravity 19(20), 5235-5247 

(2002). 

[FKNN06] S. Frittelli, C. Kozameh, E. T. Newman and P. Nurowski, Differential equations and 
Cartan connections, in Topics in mathematical physics, general relativity and cosmol- 
ogy in honor of Jerzy Plebanski, pages 193-200, World Sci. Publ., Hackensack, NJ, 
2006. 

[GN07] M. Godlihski and P. Nurowski, GL(2,R) geometry of ODE's, 2007. 

[God08] M. Godliriski, Geometry of third-order ordinary differential equations and its appli- 
cations in general relativity, PhD thesis. University of Warsaw, 2008. 

[Gra87] C. R. Graham, On Sparling's characterization of Fefferman metrics, Amer. J. Math. 
109(5), 853-874 (1987). 

[Lee86] J. M. Lee, The Fefferman metric and pseudo-Hcrmitian invariants. Trans. Amer. 

Math. Soc. 296(1), 411-429 (1986). 
[Lie05] H. Liebmann, Beriihrungstransformationen, in Enzyklopadie der Mathematischen 

Wissenschaften, Band III, Teil D7, pages 490-492, B.C. Teubner, 1905. 
[LS91] S. Lie and G. Scheffers, Geometric der Beriihrungstransformationen, B.C. Teubner, 

Leipzig, 1891. 

[NS03] P. Nurowski and G. A. Sparling, Three-dimensional Cauchy-Ricmann structures and 
second-order ordinary differential equations. Classical Quantum Gravity 20(23), 4995— 
5016 (2003). 



CAUSAL GEOMETRIES AND THIRD-ORDER ORDINARY DIFFERENTIAL EQUATIONS 23 

[Nur05] P. Nurowski, Differential equations and conformal structures, J. Gcom. Phys. 55, 
19-49 (2005). 

[01v93] P. J. Olver, Applications of Lie groups to differential equations, volume 107 of Grad- 
uate Texts in Mathematics, Springer- Verlag, New York, second edition, 1993. 

[01v95] P. J. Olver, Equivalence, invariants, and symmetry, Cambridge University Press, 
Cambridge, 1995. 

[PR88] R. Penrose and W. Rindler, Spinors and space-time. Vol. 2, Cambridge Monographs 

on Mathematical Physics, Cambridge University Press, Cambridge, second edition, 

1988, Spinor and twistor methods in space-time geometry. 
[SY98] H. Sato and A. Y. Yoshikawa, Third order ordinary differential equations and Legcndre 

connections, J. Math. Soc. Japan 50(4), 993-1013 (1998). 
[TodOO] K. P. Tod, Einstein-Weyl spaces and third-order differential equations, J. Math. Phys. 

41(8), 5572-5581 (2000). 
[Wil06] E. Wilczynski, Projective differential geometry of curves and ruled surfaces, Teubner, 

Leipzig, 1906. 

[Wiin05] K. Wiinschmann, Ubcr Beruhrungsheingundcn bei Differentialbleichungcn, PhD the- 
sis, Greifswald, 1905. 



